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ABSTRACT 

The curvature of a relativistic blast wave implies that its emission arrives to 
observers with a spread in time. This effect is believed to wash out fast vari- 
ability in the light curves of GRB afterglows. We note that the spreading effect 
is reduced if emission is anisotropic in the rest-frame of the blast wave (i.e. 
if emission is limb-brightened or limb-darkened). In particular, synchrotron 
emission is almost certainly anisotropic, and may be strongly anisotropic, 
depending on details of electron acceleration in the blast wave. Anisotropic 
afterglows can display fast and strong variability at high frequencies (above 
the 'fast-cooling' frequency). This may explain the existence of bizarre fea- 
tures in the X-ray afterglows of GRBs, such as sudden drops and flares. We 
also note that a moderate anisotropy can significantly delay the 'jet break' in 
the light curve, which makes it harder to detect. 

Key words: gamma ray: bursts; shock waves; radiation mechanisms: non- 
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1 INTRODUCTION 

GRB afterglows are likely produced by relativistic blast waves propagating from the center 
of the explosion. This model is, however, challenged by recent observations. In particular, 
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the Swift satellite revealed several puzzling features in the X-ray afterglow. It observed an 
early plateau stage and flares with fast rise and decay times (Nousek et al. 2006; Burrows 
et al. 2005). Less frequent but even more bizarre are sudden drops in the X-ray light curve 
(as steep as t"^" in GRB 070110, Troja et al. 2007). These behaviors are inconsistent with 
the standard model of afterglow production. 

Can the emission from the forward or reverse shock of the blast wave show strong vari- 
ations on timescales At ^ t? It is usually argued that this is impossible: the spherical 
curvature of the emitting surface (of radius R and Lorentz factor F) implies a spread in 
arrival times of its emission, which washes out variability on timescales shorter than 

For a relativistic blast wave, this duration is comparable to the observed time passed since the 
beginning of the explosion, r ~ t. This appears to prohibit any rapid and strong variations 
in the light curve (see loka et al. 2005 for discussion). 

Therefore, the observed fast variability in afterglows is usually associated with additional 
emission from radii much smaller than the blast-wave radius. This model invokes a late 
activity of the central engine (Zhang et al. 2006). The material ejected at large t and emitting 
at radii R <^ F^tc will have t <^t and can produce flares with At ^ t. Note however that 
(i) it is unclear in this model why the observed flares have the approximately universal 
At/t ~ 0.1 (Chincarini et al. 2007; Lazzati & Perna 2007), (ii) the very steep drops at the 
end of some plateaus can hardly be explained by this model unless it assumes that the entire 
plateau is produced at small radii inside the ejecta and the emission from the blast wave is 
negligible (Kumar, Narayan & Johnson 2008). 

Another difficulty for GRB theory is that many afterglows lack the predicted 'jet breaks' 
(Burrows & Racusin 2006; Sato et al. 2007): only a small fraction of afterglow light curves 
show a clear achromatic break that is expected from jets (Willingale et al. 2007), l| Some 
bursts show X-ray light curves extending for tens to hundreds of days with a constant 
temporal slope (Grupe et al. 2007). The interpretation of these observations is difficult and 
often leads one to assume large jet opening angles, implying in some cases extremely high 
energy for the explosion (Shady et al. 2007). 

An implicit assumption in the general discussion of these puzzling features is that the 
emission is isotropic in the rest frame of the relativistically moving source (however, see 



^ Many afterglows show chromatic breaks, which occur either in the X-ray or in the optical, but not in both bands. 
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Lyutikov 2006). In this paper, we discuss the effects of a possible anisotropy and suggest 
that they can help explain observations. In Section[2]we write down a general formula for the 
observed flux from a flashing sphere when the emission is anisotropic in the source rest frame. 
In Section [3] we list the consequences of anisotropy for the curvature effect, the jet break 
in the afterglow light curve and the size of the radio image of the blast wave. In Section H] 
we consider the standard radiative mechanism of afterglows - synchrotron emission - and 
discuss its anisotropy. The results are summarized in Section |5l 



2 LIGHTCURVE FROM A FLASHING SPHERE 

Let an energy E (measured in the lab frame) be instantaneously emitted by a sphere of 
radius R, which is expanding with a Lorentz factor F and a velocity P = v/c = (1 — F^^)^^^. 
Distant observers will see the emitted radiation extend over a range of arrival times t due 
to the curvature of the emitting surface. The observed light curve L{t) from the flashing 
sphere can be thought of as a Green function of afterglow emission or 'response function' 
that describes the curvature effect. It depends on the intrinsic angular distribution of the 
source intensity. Let 6 be the photon angle with respect to the radial direction in the local 
rest frame of the emitting sphere. The intrinsic (comoving) angular distribution of emission 
(per unit solid angle) can be described by function A{6) normalized by 

J A{e) dQ = An. (2) 

Isotropic emission in the comoving frame corresponds to A{9) = 1. The photon angle in the 
static lab frame 6 is related to 6 by 
cose + P 

cos 9 = :; . (3) 

A distant observer will first see photons emitted along the line of sight with 6 = Q = 0. 
Let to be the arrival time of these first photons. Photons received at a later time t come 
from a larger co-latitude B on the emitting sphere, related to t by 
R 

t-to = -(1 -cosG). (4) 

A time interval 6t corresponds to a ring 5cosG = c6t/R on the sphere. The true energy 
emitted by this ring is 6E = 6 cosQ. The ring is viewed at angle with respect to 
its normal, and the observed photons are Doppler-boosted in energy by the factor V = 
F~^(l — /3cos6)~^. Doppler effect also compresses the solid angle of emission by the factor 
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of V^^. Together with anisotropy A{6), Doppler effect determines the amphfication factor 
for the apparent isotropic equivalent of emitted energy, 

6E^,, = A{e)V'^-^. (5) 
The observed luminosity is L = 6Eg^pp/6t, which yields 

L = ^^^^ (6) 

2i?r4(i -/3cose)3' ^ ^ 

where 9 and 9 are related by Equation ([3]). Substitution of cos6(t) = 1 — c(t — to) /R and 
6'[0(t)] gives an explicit expression for the light curve produced by the flashing sphere. These 
equations simplify in the limit F ^ 1, 
r - (t - to) 

cos^ = — — — — , (7 
r + [t- to) 

L^?lm{^.i^)-\ (s, 

where r is defined in Equation ([T]). 



3 SOME CONSEQUENCES OF ANISOTROPY 

An isotropic source, A{6) = 1, after Doppler transformation to the static frame emits 75 per 
cent of the energy within Obeam = l/f - Let us now consider an anisotropic source, A{6) ^ 1. 
We will assume that the anisotropy has a front-back symmetry, A{Q) = Aln — 6); then the 
net momentum of emitted photons vanishes in the source frame. 

Consider, for instance, 'limb-darkened' emission, which is weak near 6 = 7i/2 and strong 
near 6 = 0,7t. Doppler transformation to the lab frame strongly amplifies the radiation with 
^ ^ and weakens radiation with 6* ^ vr. As a result, a bright narrow beam is created, so 
that 75 per cent of energy is now concentrated within Gbeam = (^r)^^. Here A; > 1 is a 
measure of the enhanced beaming of radiation in the lab frame. The beam < 0bcam is 
emitted with 9 < ^^bcam in the source frame, and one can show that 6'beam is related to k by 

/ l + COS^beam ^ 
V 1 - COS 6'beam 

The increased beaming in the lab frame {k > 1) due to limb-darkening in the source frame 
(cos 6'beam > 0) has scvcral observational consequences that we list below. 



3.1 Curvature effect 

The curvature effect is expected to control the observed hght curve if the source power 
suddenly drops. The observed luminosity L{t) responds to the drop with a delay according 
to Equation ([8]). If the emission is isotropic in the source frame, the delay timescale is r ~ t, 
and the steepest possible decay is L(t) oc t'" with a = 3 (e.g. Kumar & Panaitescu 2000). 
Since limb-darkening of the source implies stronger beaming, Bbeam = (^r)~^ instead of F"^, 
most of the energy is radiated on a shorter timescale, r/k"^, and the slope a is much steeper. 

We will discuss this effect in more detail in Section 4. It turns out that a similar conclusion 
holds for the opposite, limb-brightened, type of anisotropy, when emission is strong near 
9 = n/2 and weak near 9 = 0,n. 

3.2 Jet break 

GRB jets are likely to have a small opening angle 0jct <^ 1, which reduces their energy 
Ejet ~ (0fet/2)-E'iso (here i^iso is the isotropic equivalent of the jet kinetic energy). A break 
should be observed in the afterglow light curve at moment tjet when the relativistic beaming 
angle 6beam becomes larger than ©jet (Rhoads 1997), and the value of Gjet may be inferred 
from observed tjet- If the afterglow source is limb-darkened so that 6beam = {kT)~^, the 
jet-break condition becomes T ^ (A;6jet)~^, i-e. effectively 6jet is replaced by kQ^ct- The 
standard light-curve analysis can only give the value of fcGjet, which overestimates the true 
Gjet by the factor of k. The true Ej^t for a limb-darkened jet is reduced by the factor of k~^ 
compared with the usual estimate. 

Limb-darkening also implies a significant delay in tjet. For example, consider a blast wave 
decelerating in a uniform medium. Its Lorentz factor decreases as F oc t^^^^. The jet break 
occurs when F (fcBjet)"^, and this moment is delayed by the factor of k^^^. The usual 
expression for tjet then becomes 

/TP X 1/3 X 8/3 

*M«*"'(%i) Of) d. (10) 

Similarly, for a blast wave decelerating in a wind medium, F oc t~^^^ and hence tjet oc fc^. Even 
a moderate anisotropy (e.g. k = y/S, which corresponds to limb- darkening with 6'beam ~ 60° 
in the source frame) can delay the jet break by a large factor (~ 4.3 in a uniform medium, 
~ 9 in a wind). This could be enough to not detect the jet break with current observational 
capabilities as the afterglow is dim at late times and the spectral coverage is incomplete to 
test achromaticity. 
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3.3 Apparent size of the radio afterglow source 

VLBI observations provided the angular size of the radio image of a few GRB afterglows 
(Frail et al. 1997; Taylor et al. 2005) which helps constrain the ratio of the blast wave energy 
to the density of the environment. The best data have been obtained for GRB 030329 and 
seem to favor a blast wave in a uniform medium with Ejct/n ~ (1 — 5) x 10^^ erg cm^ and 
Gjet ~ 0.1 rad (Pihlstrdm et al. 2007). 

The apparent size of the afterglow source is given by (e.g. Oren, Nakar & Piran 2004) 

R^C^RQ, e = min(ebeam;ejet), (11) 

where R is the radius of the emitting shell. With increased beaming due to limb-darkening, 
©beam = ^ud R_i IS reduccd by the factor k~^. For a blast wave in a uniform medium, 

a derivation similar to that in Oren et al. (2004) gives the relation between the ratio of the 
true jet energy to the external density and observed R±, t and tjet, 

^jet , .4 ( R± \ ' \ €t t''"' (before jet break) 3 

^ 10 k X erg cm"^ (12) 

n KQxW^cmJ y ^-3 (after jet break) 

One can show that k drops out from the similar relation derived for blast waves in wind 
media, i.e. in that case limb-darkening does not affect the relation. 



4 ANISOTROPY OF SYNCHROTRON EMISSION 

Afterglow is commonly interpreted as synchrotron emission. Its anisotropy naturally results 
from a preferred orientation of the magnetic field B. Magnetic fields inside GRB jets are 
generally expected to be transverse to the jet direction, as radial expansion quickly suppresses 
the longitudinal component. Internal or external shocks can generate magnetic fields only 
in the shock plane (e.g. Medvedev & Loeb 1999). Thus, in various models of afterglow 
productioij^ it is reasonable to suppose that the magnetic field in the source is perpendicular 
to its velocity, B = B_l. In addition, we assume that B_l is tangled on a scale much smaller 
than i?/r, so that a distant observer will see a superposition of emissions from many domains 
with random orientations of B^. This assumption is motivated by the low polarization in 
observed afterglows, typically less than a few per cent (e.g. Covino et al. 1999). 



^ At present, the origin of afterglow emission is unclear. The standard forward-shock model is in conflict with data, and it is 
possible that the afterglow is produced by a long-lived reverse shock (Uhm & Beloborodov, 2007; Genet, Daigne & Mochkovitch, 
2007). 
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The anisotropy of synchrotron emission may be expected to be moderate if the emitting 
electrons have an isotropic distribution (see e.g. calculations by Granot et al. (1999) for three 
possible geometries of the magnetic field). Even in this case, anisotropy is present because B 
is confined to a plane. After averaging over random directions of B^, one finds the angular 
distribution of emitted power per electron, 

^ = ^A,ie), with Aoie) = Ul + cos'0). (13) 

Here Pq = (crTc/67r)7^P^ is the power of synchrotron emission per electron and is Thom- 
son cross section. The resulting radiation is limb-darkened. 

In reality, the electron distribution may not be isotropic: electrons may be preferentially 
accelerated along or perpendicular to the magnetic field, depending on the acceleration 
mechanism. For instance, the details of electron acceleration in relativistic shocks remain 
uncertain, despite significant progress in numerical simulations (e.g. Hededal et al. 2004; 
Spitkovski 2008; Nishikawa et al. 2009), and other mechanisms are possible. We therefore 
consider both types of electron anisotropy. 

Let a be the pitch-angle of an electron with respect to the magnetic field. We will 

describe the distribution of electron directions fie by the function /(a), normalized by 

/ /(«) dQf, = Att. Synchrotron emission from each electron is strongly beamed along its 

velocity, and together the electrons emit radiation with angular distribution, 

dP 3Po . 2 rr \ r-tA\ 
^ = — sm «/(«). (14) 

The average power per electron is now given by, 

P = r/Po, V = -^ [ sin' a f{a) d^,. (15) 

Equation (JT^ describes the angular distribution relative to the local magnetic field. The 
angle a between the magnetic field and the observer's line of sight is given by cos a = 
sin ^ cos where < < 27r depends on the orientation of B = B^. Using Equation (IHj) 
and averaging over random directions of B^, one finds that the synchrotron emission has 
the following angular distribution, 

dP nPn 3 1 / \ 

— = ^A{e), with A{e) = —— ( 1 - sin' cos' 0) / [arccos(sin^cos0)] (16) 
aiZ 47r zf] ZTT Jo ^ ' 

Let us consider two toy models: 
j\{oL) oc [g? + sin' and /2(tt) oc {c? + cos'aj , (17) 
where a defines a characteristic beaming angle of the electron distribution. These expressions 



represent two opposite cases where electrons are preferentially accelerated along B (distri- 
bution /i) and perpendicular to B (distribution /2). Both distributions becomes isotropic 
if a ^ 1. The opposite limit a ^ 1 describes the maximum possible anisotropy: /i = 6 {a) 
and /2 = S{a — 7i/2). By varying the parameter a from oo to 0, one can explore the effect 
of increasing electron anisotropy on synchrotron emission. 

The angular distributions Ai{9) and ^2(6*) produced by /i and /2 are shown in Figures [1] 
and [21 The plot of Ai{6) resembles a butterfly with a half-angle ~ a around 9 = 7r/2, i.e. the 
source is limb-brightened. Therefore, the angular distribution of radiation in the fixed lab 
frame has a sharp peak at B = F^^ when a ^ 1. By contrast, A2{6) is concentrated near 
6 = 0,TT with a half-angle ~ a; the resulting limb-darkening remains, however, finite even if 
a — )■ 0. In this limit, one finds 

^^^'^ = dre (1^) 

The source is weakened at 6 = 71 /2 by the modest factor of 2/7r compared with isotropic 
emission^ Note also the difference in the emission power P for the two distributions. For 
a ^ 1 one obtains 
3 3 

r/l^-a^ ^2~-. (19) 

The anisotropy of emission A{6) ^ 1 can have a strong impact on the afterglow light 
curves as discussed in section 3. In particular, the curvature effect, which is described by the 
light curve from a flashing sphere (the response function. Equation ([8])), is changed. Figure [3] 
shows the response function for A\{Q^ and ^2(6*) with a = 0.03. Compared to the isotropic 
case, A{6) = 1, the emitted pulse becomes very narrow if A{6) = Ai{6), i.e. for the model 
where electrons are preferentially accelerated along the magnetic field. Photon arrival times 
then concentrate near a particular t — to ^ t (which corresponds to a particular = F^^) 
because the limb-brightened radiation is mainly emitted near 6* = 7r/2. In the model with 
angular distribution ^2(6'), the profile of the response function is steeper than in the isotropic 
case, but can never be as narrow as for Ai{6), even in the limit of a 0. 

As a simple illustration, consider a spherical thin shell with constant emission power Eq, 
which is moving with Fq = 300, and suppose that its emission suddenly cuts off at radius 
Rcut = 6 X 10^® cm. Figure m shows the produced bolometric light curve. It depends on the 



^ Electron distribution with /2 = S{a~ n/2) is a disk in momentum space, with the axis along B. After averaging over random 
orientations of the disk axis in the transverse plane (random directions of B) a strong anisotropy is found in A2{9), with a 
sharp peak in the longitudinal directions 9 = 0, tt. However, a significant 'wing' of emission remains present at large 8 ~ it/2. 
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Figure 1. Diagram of angular distribution of radiation Ai{6) (electrons accelerated preferentially along B) measured in the 
source rest frame (left) and transformed to the observer frame using F = 10 (right). Solid curves correspond to a = 0.1 and 
dashed to a = 0.03. Long-dashed lines in the right panel show the cone of opening angle l/F. Rotation of the shown curve 
about the horizontal axis gives the 3-dimensional diagram. 
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Figure 2. Same as in Figure [T] but for the angular distribution A2{0) (electrons accelerated preferentially perpendicular to B). 

intrinsic anisotropy of the source, ^4(6'). We show three cases: isotropic emission ^4(6') = 1, 
Ai{6) and A2{6) (same as in Fig. 3). The hmb- brightened emission ^1(6*) produces a very 
steep decay in the hght curve. 

We conclude that extremely fast variations in the light curve, e.g. short flares or steep 
drops, may be observed in synchrotron afterglows if electrons are preferentially accelerated 
along B, as the response function can be arbitrarily narrow for a ^ 1. Note that the 
emission is limb-brightened in this case, i.e. the situation is opposite to what was considered 
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(t-g/T 

Figure 3. Lightcurve from a flashing sphere with three difl'erent angular distributions of emission in the plasma rest frame: 
isotropic (dotted curve), Ai{B) with a = 0.03 (solid curve), and A2{0) with a = 0.03 (dashed curve). The units for t and L are 
indicated on the axes (r = R/2r^c); with these units the area under each curve is unity. The first photon from the flashing 
sphere reaches the observer at t = tg. The dashed curve has the maximum of about 13 at t = to- 
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Figure 4. Bolometric light curve from a shell whose emission suddenly cuts off (see text). 

in section 3, and the description using k > 1 does not apply. Limb-brightening has httle or 
no effect on tjet or R_i, in contrast to the hmb-darkened model of section 3. 

Another implication of the preferential electron acceleration along B is the reduction 
of synchrotron emissivity by a factor of ~ (cf. r/i in Equation f|T9|) ). All synchrotron- 
emission formulae contain only the component of B perpendicular to the electron velocity, 
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which equals ~ aB for a <^ 1. Then the effective that would be inferred from the data 
using isotropic models will underestimate the real by a factor ~ a~^. 

The synchrotron model with angular distribution A2{9), does not predict significant 
changes in the afterglow light curve, because limb-darkening is never strong for synchrotron 
emission, regardless of a. A moderate change in tjet and a less pronounced jet break may be 
expected compared with the case of isotropic emission. 

5 DISCUSSION 

The usual assumption of isotropic emission in the rest frame of the blast wave is likely 
to be invalid. Even the standard synchrotron model with isotropic electron distribution 
produces anisotropic, limb-darkened radiation (Section H]). This fact is a consequence of the 
preferential orientation of the magnetic field in the blast wave. Strong limb-brightening is 
also possible if the radiating electrons are preferentially accelerated along the magnetic field. 

Anisotropy may resolve a few puzzles encountered in afterglow modeling: 
(i) The usual argument that the curvature effect filters out fast variability, prohibiting strong 
variations in the light curve on timescales At < r = i?/2r^c, is not valid for anisotropic 
emission. An anisotropic variable spherical source can produce fast changes in the light 
curve, similar to observed bizarre features in GRB afterglows. This result holds for both 
limb-darkened and limb-brightened types of anisotropy. It suggests that the X-ray flares 
observed by Swift with At/t<0.1 do not necessarily imply an additional component of 
internal origin. Instead they may be produced, e.g. by the reverse shock in the blast wave, 
whose emission may suddenly brighten and weaken as the reverse shock propagates into the 
inhomogeneous ejecta of the explosion. This model may also explain sudden steep drops in 
the afterglow light curve as observed in GRB 070110 (Troja et al. 2008). This explanation 
assumes that the X-ray radiating particles are cooling fast compared with the jet expansion 
timescale, as slow cooling would suppress short time-scale variations of the source luminosity. 

Examples of such short times-cale behaviors are given by the toy model in Figure [51 It 
shows the synchrotron emission produced by a thin shell with Lorentz factor T{R) = Tq = 
300 at i? < i?doc = 3 X 10^'' cm and T{R) = Tq {R/ Rdec)"^^'^ at i? > Rdcc- This approximately 
describes a blast wave decelerating in a uniform medium. The shell is assumed to radiate 
with bolometric power proportional to the dissipation rate in the blast wave, which gives 
E{R) = Eo {R/Rdecf at R< Rdec and E{R) = Eq (R/Rdec)'^ &t R > R^ec- A realistic blast 
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Figure 5. Bolometric light curve for a toy afterglow model (see text). The result is plotted for three cases: isotropic emission 
in the rest frame of the blast wave A{0) = 1 (solid curve), limb-brightened synchrotron emission Ai{9) with a = 0.03 (dotted 
curve; sec Section 4 for the description of the synchrotron model) and limb-darkened synchrotron emission A2{9) with a = 0.03 
(dashed curve). 



wave has two shocks - forward and reverse - and both can produce a long-hved afterglow. Our 
toy model may describe the emission from either shock, although it is very much simplified. 
To illustrate the curvature effect on variabihty, we add two features: a sudden brief increase 

in E{R) aX R = 3-Rdcc (which simulates a flare) and the abrupt cutoff of E{R) at 15-Rdcc- 
For comparison, we show the light curves produced for three cases: isotropic emission in the 
rest frame of the shell, limb-brightened emission and limb-darkened emission described in 
Section 4. 

(ii) If a relativistic source is limb-darkened, most of its emission in the fixed lab frame 
is confined within an angle smaller than F"^. This effect suggests a possible explanation 
for the lack of jet-break detections in GRBs, as the increased beaming significantly delays 
the jet break in the observed light curve (Section 3.2). We also discussed in Section 3.3 the 
consequences of such anisotropy for the apparent size of the radio afterglow source. Although 
the strong limb-darkening appears to be impossible for standard synchrotron afterglows, it 
may be possible for a different radiative mechanism. For example, limb-darkening may be 
expected for the jitter mechanism (Medvedev & Loeb 1999), as the electrons are preferen- 
tially accelerated perpendicular to the shock plane and radiate preferentially in the radial 
direction. 

While this paper was focused on afterglow, the source of prompt GRB emission may also 
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be intrinsically anisotropic. This may impact models that propose the curvature effect to 
control the steep X-ray decay at the end of the prompt emission (see e.g. Genet & Granot 
2009; Zhang et al. 2009). The effect can be seen in Figure 4. Anisotropy of the prompt 
emission may also change the optical depth of the source to high-energy photons, r^^, as 
the cross section for 77 reaction strongly depends on the angle between photons. This may 
affect the constraints on the Lorentz factor of the jet that are inferred from r^^ < 1. The 
effect is especially strong for emission without front-back symmetry in the source frame; 
such asymmetric emission would be a more radical assumption compared with the ordinary 
limb-brightening or limb-darkening considered in this paper. 



ACKNOWLEDGMENTS 

We thank Jonathan Granot and the referee for comments on the manuscript. AMB was 
supported by NASA Swift grant, Cycle 4. ZLU is supported by WCU program (R32-2009- 
000-10130-0) of NRF/MEST of Korea. FD and RM were supported by the French space 
agency (CNES). 



REFERENCES 

Burrows D.N., Racusin J., 2006, NCimB, 121, 1273 
Burrows D.N. et al., 2005, Sci, 309, 1833 
Chincarini G. et al., 2007, ApJ, 671, 1903 
Covino S. ct al., 1999, A&A, 348, LI 
Frail D.A., 1997, Nature, 389, 261 

Genet F., Daigne F., Mochkovitch R., 2007, MNRAS, 381, 732 

Genet F. & Granot J., 2009, MNRAS, 399, 1328 

Granot J., Piran, T., Sari, R., 1999, ApJ, 527, 236 

Grupe D. et al., 2007, ApJ, 662, 443 

Hededal, C.B. et al., 2004, ApJ 617, L107 

loka K., Kobayashi S., Zhang B., 2005, ApJ, 631, 429 

Kumar R, Narayan R., Johnson J.L., 2008, MNRAS, 388, 1729 

Kumar P., Panaitescu, A., 2000, ApJ, 541, L51 

Lazzati D., Perna R., 2007, MNRAS, 375, L46 

Lyutikov, M., 2006, MNRAS, 369, L5 

Nishikawa, K.-I. et al., 2009, ApJ 698, LIO 

Nousek J.A. et al., 2006, ApJ, 642, 389 

Oren Y., Nakar E., Piran T., 2004, MNRAS, 353, L35 

Pihlstrom Y.M., Taylor G.B., Granot J., Doeleman S., 2007, ApJ, 664, 411 
Rhoads J., 1997, ApJ 487, LI 
Sato G. et al., 2007, ApJ, 657, 359 



14 



Shady P. ct al., 2007, MNRAS, 380, 1041 

Spitkovski, A., 2008, ApJ 682, L5 

Taylor G.B. et al., 2005, ApJ, 622, 986 

Troja E. et al., 2007, ApJ, 665, 599 

Uhm Z.L., Beloborodov A.M., 2007, ApJ, 665, L93 

Willingale R. et al., 2007, ApJ, 662, 1093 

Zhang B. et al., 2006, ApJ, 642, 354 

Zhang B.-B. et al., 2009, ApJ, 690, LIO 



